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Abstract —In two-tiered Wireless Sensor Networks (WSNs) 
relay node placement Is one of the key factors impacting the 
network energy consnmption and the system overhead. In this 
paper, a novel connectivity-aware approximation algorithm for 
relay node placement in WSNs is proposed to offer a major step 
forward in saving system overhead. Specifically, a uniqne Local 
Search Approximation Algorithm (LSAA) is introduced to solve 
the Relay Node Single Cover (RNSC) problem. In this proposed 
LSAA approach, the sensor nodes are allocated into groups and 
then a local Set Cover (SC) for each group is achieved by a 
local search algorithm. The union set of all local SCs constitutes 
a SC of the RNSC problem. The approximation ratio and the 
time complexity of the LSAA are analyzed by rigorous proof. 
Additionally, the LSAA approach has been extended to solve 
the relay node double cover problem. Then, a Relay Location 
Selection Algorithm (RLSA) is proposed to ntlllze the resulting 
SC from LSAA in combining RLSA with the minimum spanning 
tree heuristic to build the high-tier connectivity. As the RLSA 
searches for a nearest location to the sink node for each relay 
node, the high-tier network bnilt by RLSA becomes denser than 
that by existing works. As a result, the number of added relay 
nodes for building the connectivity of the high-tier WSN can be 
significantly saved. Simulation results clearly demonstrate that 
the proposed LSAA outperforms the approaches reported in 
literature and the RLSA-based algorithm can noticeably save 
relay nodes newly deployed for the high-tier connectivity. 

Index Terms —wireless sensor networks, relay node placement, 
geometric disc cover, network connectivity, approximation algo¬ 
rithm. 

I. Introduction 

IRELESS Sensor Networks (WSNs) have been one of 
the most important research areas since last decade 
due to its tremendous application potentials in civilian and 
defense related applications such as battlefield surveillance, 
environmental monitoring, biomedical observation, industrial 
automation and other fields Q-ISl. WSNs typically composed 
of many low-cost and low-power homogenous or heteroge¬ 
neous sensor nodes, which can perform sensing, simple com¬ 
putations, and short-range wireless communications. In many 
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WSN applications, the lifetime of network is limited due to the 
constraints in energy resources and accessibility of the actual 
sensor nodes a-a. The two-tiered network architecture has 
been proposed to extend the network lifetime in WSNs. In 
this approach a small number of nodes with ample power and 
suitable wireless communication radius are placed to serve as 
the relay nodes. The relay nodes function like the cluster heads 
to collect sensed information from one-hop-neighbor sensor 
nodes and transmit the data to the sink node a-il. This is 
an energy-efficient approach since it reduces the energy and 
resource consumption of each sensor node and relies on the 
mesh networks of relay nodes to transmit data from sensor 
nodes. This highlights the key role of relay nodes and the 
importance of investigating their placement considering the 
required energy consumption and other overheads. 

The research problems on relay node placement in two- 
tiered WSNs can be divided into two subgroups, 1) the Relay 
Node Cover (RNC), which is mathematically equivalent to the 
Geometric Disc Covering (GDC) problem, and 2) the network 
connectivity (NC). Different published works have investigated 
the problems of GDC and NC independently, however, this 
paper considers the two subproblems jointly, i.e., the relay 
nodes for covering the low-tier WSN are optimally deployed 
to facilitate the connectivity of the high-tier WSN. 

This paper first focuses on the GDC problem and proposes a 
novel Local Search Approximation Algorithm (LSAA) to the 
Relay Node Single Cover (RNSC) problem. Different from 
noi-iEi], LSAA is a two-phase algorithm, in which the first 
phase allocates the sensor nodes into groups and the second 
phase searches a Set Cover (SC) for the GDC problem. The 
global SC consists of all local SCs that are found by a new 
local search algorithm called-Neighbor First Weighted Greedy 
Algorithm (NFWGA), with each SC covering one group. 
Here, the approximation ratio and the time complexity of the 
LSAA are studied rigorously. Additionally, a unique LSAA 
is proposed and evaluated to address the relay node Double 
Cover problem. 

This paper then studies the NC problem that is NP-hard 
and proposes a Relay Location Selection Algorithm 
(RLSA)-based NC algorithm. In this approach, RLSA first 
selects an optimal deployed location (nearest to the sink node 
in Euclidean distance) for each relay node, by utilizing the SC 
returned by LSAA, in building the high-tier network connec¬ 
tivity based on the Minimum Spanning Tree (MST) heuristic 
ED. As a result of the RLSA outstanding performance, the 
high-tier network becomes denser and thus fewer relay nodes 
are added to the high-tier network for connectivity. 

In summary, based on results and comparisons with the 
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published results, this work is dehning a new and important 
direction in relay node placement for two-tiered WSNs. Some 
of the signihcant contributions of this new approach are listed 
below: 

• Proposing a novel solution for the GDC problem, called 
Local Search Approximation Algorithm, which consists 
of two steps of hrst allocating the sensor nodes into in¬ 
dependent groups and then proposing NFWGA to search 
a local SC for each group. Detailed and rigorous analysis 
on the approximation ratio and the time complexity of 
the LSAA are also provided in this paper. 

• Proposing RLSA that selects one optimal deployed loca¬ 
tion for each relay node returned by LSAA. As a result, 
the RLSA-based NC algorithm saves newly added relay 
nodes for maintaining the network connectivity. 

• Extensive simulations are conducted to demonstrate the 
efficiency of the algorithms proposed in this paper. 

This paper is organized as follows. Section 11 reviews related 
work. Section III details the problem formulation. Section IV 
presents the approximation algorithm for the GDC problem. 
Section V presents the RLSA-based NC algorithm. Simulation 
results are presented in Section VI to demonstrate the superior 
performance of this work. Section VII hnally concludes this 
paper. 

11. Related Work 

A. GDC Problem 

In two-tiered WSNs, the placement of relay nodes is a sig¬ 
nihcant factor in energy consumption and networks overhead. 
In order to address the NP-complete nature of the relay node 
placement problem ifTOl . existing research works mainly focus 
on three approaches of the Shift Strategy, the Grid Strategy 
and the Set-Covering Strategy. In this section, let n, I and e 
denote the number of sensor nodes, the shift parameter dehned 
by the Shift Strategy and any positive constant, respectively. 

The Shift Strategy: D. S. Hochbaum ifTTl ifTSll studied 
the GDC problem and proposed a Polynomial Time Approx¬ 
imation Scheme (PTAS) based on the Shift Strategy, which 
divides the deployed held into stripes with the same width 
and partitions these stripes into different groups by shifting 
the stripes, then solves the GDC problem in each stripe. 
The approximation ratio and time complexity of their scheme 
are bounded by (l+l/l)'^ and 0(P [v/2/] with 

Z > 1. T. Feder and D. Greene Ha, independently, Gonzalez 
na considered the related problem and proposed a (1 -f 1/Z)- 
approximation algorithm, which solves the problem by divid¬ 
ing the deployed held into stripes, with I > 1, and the time 
complexity of this algorithm is J. Tang 

et ai. na proposed two approximation algorithms based on 
the PTAS proposed in im ca to solve the relay node single 
cover problem and relay node double cover problem based on 
the Shift Strategy, with the approximation ratio of 4 and 9/4, 
respectively. A. Srinivas et al. ESI considered the problem 
of the maintenance of a mobile backbone using the minimum 
number of backbone nodes, and they also proposed algorithms 
for the GDC problem based on the PTAS proposed in EH 
Ell, where the relay node had no mobility. However, the time 


complexity of the Shift Strategy grows exponentially with the 
shift parameter and the strategy is rather time-consuming when 
the shift parameter is large. 

The Grid Strategy: M. Franceschetti et al. EtI Esl 
proposed a grid strategy, which divides the deployed held 
into square meshes and the relay nodes can only be placed 
at the vertices of square meshes. They hrst analyzed the ap¬ 
proximation factor of grid strategy under different disc radius. 
Then they combined the Grid Strategy with the Shift Strategy, 
which leads to a family of algorithm with a performance ratio 
of cr(l -f 1/Z)^, with cr G {3,4, 5,6} and I > 1. The poor 
approximation ratio of the Grid Strategy stems from the fact 
that all placement locations in one grid are hxed. 

The Set-Covering Strategy: In the Set-Covering Strategy, 
each relay node is represented by the sensor nodes covered 
by it, correspondingly, the Set-Covering Strategy searches a 
minimum set cover of the given sensor nodes. Bronnimamm 
and Goodrich ED studied the set covering problem of dual 
VC-dimension and proposed an algorithm with a constant 
approximation ratio, which is not specihed yet. Then they 
gave several applications of their method to computational 
geometry including the GDC problem. K. N. Xu et al. Il20l hrst 
modeled GDC problem with a minimum set covering problem 
in WSNs, and proposed a recursive algorithm. K. Ah et al. 
m considered the related problem and proposed a weighted 
greedy algorithm, with an approximation ratio of ln(n) -f 1. 
Obviously, the Set-Covering Strategy yields an approximation 
ratio which grows logarithmically with the scale of the network 
and suffers from the scalability issue. 

Unlike Set-Covering and Grid strategies that solve the whole 
GDC problem at a time, LSAA hrst decomposes the whole 
GDC problem into several GDC subproblems with smaller 
size (which is also different from the Shift Strategy that simply 
divides the deployed held into stripes), and then combines the 
solutions to each GDC subproblem to form the solution to 
the whole GDC problem. By the novel decomposition-based 
relay deployment method, LSAA has a small approximation 
ratio and a low time complexity. 

B. Network Connectivity Problem 

Extensive works have been done ED 123-133 for the 
network connectivity problem and can be classihed into two 
branches, i.e., single-tiered network connectivity and two- 
tiered network connectivity. 

In the single-tiered network, sensor nodes act as both 
transmitter and relay node. G. Lin and G. Xue 123 modeled 
the relay node placement problem in the single-tiered network 
with the Steiner Minimum Tree with Minimum number of 
Steiner Points and bounded edge length (SMT-MSP) problem. 
The SMT-MSP problem was proved NP-hard and then solved 
by a 5-approximation algorithm. Chen et al. 123 demon¬ 
strated that the algorithm proposed in 123 was actually a 
4-approximation algorithm, and proposed a 3-approximation 
algorithm with time complexity 0(n'^), based on the so-called 
4-star. Cheng et al. 1241 also studied the same problem of 
123 and presented an improved 3-approximation algorithm 
with time complexity O(n^), based on the so-called 3-star. 
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Misra et al. im m studied the constrained relay node 
placement problem in the single-tiered network and proposed 
a polynomial time 0 (l)-approximation algorithm. 

In the two-tiered network, sensor nodes transmit the sensed 
packets to neighbor relay nodes and relay nodes forward 
the received packets to the sink node. Besides the GDC 
approximation algorithms, J. Tang et al. na also proposed 
two algorithms with approximation ratio of 8 and 4.5 for the 
connected relay node single cover problem and two algorithms 
with approximation ratio of 6 and 4.5 for the 2-connected relay 
node double cover problem, respectively. E. Lloyd and G. Xue 
first proposed a 7-approximation algorithm for the relay 
node placement problem in single-tiered network, and then 
presented a (5H-e)-approximation algorithm for the relay node 
placement problem in two-tiered network. Q. Wang et al. Il28l 
divided the relay node placement problem into two phases 
which correspond to the two problems in this paper, but mainly 
designed three heuristic algorithms to the network connectivity 
problem. D. J. Yang et al. Il29l considered the constrained 
relay node placement problem in two-tiered network and 
presented two algorithms whose approximation ratios are 0 ( 1 ) 
and 0 (ln(n)) for the connected single cover problem and 2 - 
connected double cover problem, respectively. In Il30|-|l33l, the 
relay node placement problem with the consideration of fault 
tolerance was studied. The main objective of these works was 
not only the high-tier network connectivity, but also the fault 
tolerance of the network (i.e., building at least two paths for 
each pair of sensor nodes). 

The works in GDC have been one-sided that each work such 
as iflSl ll22l - ll^ solve both GDC and NC problem separately. 
The critical flaw in these works is that they ignore the 
important relation between these two problems. In this work 
here, the novelty is that these two problems are considered 
jointly which will result in facilitating the work for the NC 
problem. 


III. Problem Formulation 

Let X = {xi,X 2 , •. •, Xn} denote a position set of n sensor 
nodes, where Xi {1 < i < n) denotes the location of the 
ith sensor node. Let Y = {yi, j/ 2 , • ■ •, 2/m} denote a position 
set of m relay nodes, where yi (1 < i < m) denotes the 
location of the /th relay node. Let a positive real number r 
denote the communication radius of sensor nodes, and R the 
communication radius of relay nodes. This paper assumes R > 
2r. 

Let p and q be two points in the plane. \\p — y|| represents 
the Euclidean distance between p and q. \f \\xi — yj\\ < r, 
1 < i < n and 1 < j < m, then the sensor node Xi can 
communicate with the relay node yj, i.e., the sensor node Xi is 
covered by the relay node yj. If jjxi — Xjj| < 2r, 1 < i,j < n, 
sensor nodes Xi and Xj are neighbors. 

The RNC problem can be further divided into the Relay 
Node Single Cover (RNSC) problem and the Relay Node 
Double Cover (RNDC) problem. 

Definition 3.1 [Relay Node Single Cover (RNSC) prob¬ 
lem]. Given a set of sensor nodes X = {xi,X 2 , ■ ■ ■ ,Xn}, 
the RNSC problem searches a set of relay nodes Y = 




Fig. 1. Possible positions for the set of sensor nodes X = 
{xi,X 2 ,x^,X 4 ,x^} are implied by the shades in different colors. The set of 
possible positions is .F = {Pi, P 2 , P3, P4, P5, Pel- 



Fig. 2. Formulating an example RNSC problem as a set covering problem. 


{ 2 / 1 , 2 / 2 , ■ • ■, 2 /m}, such that each sensor node Xi (1 < i < n) 
is covered by at least one relay node. 

To achieve a better network robustness, the RNDC problem 
is also considered in this paper and defined as follows. 

Definition 3.2 [Relay Node Double Cover (RNDC) prob¬ 
lem]. Given a set of sensor nodes X = {xi,X 2 , ■ ■ ■ ,Xn}, 
the RNDC problem searches a set of relay nodes Y = 
{ 2 / 1 , 2 / 2 , ■ ■ ■, 2 /m}, such that each sensor node Xi (1 < i < n) 
is covered by at least two relay nodes. 

The communication range of a sensor node Xi can be 
described as a disk centered at location Xi with the radius r. 
In the RNSC problem, the potential locations of relay nodes 
should be selected in the first step. For a position set of sensor 
nodes V (1 < jXl < n), if there exists a point set P such 
that ||i — pII < r, \/p G P and x G X, then P is considered 
as a possible position to deploy relay nodes. For each sensor 
node having no neighbors, its communication range could be 
selected as a possible relay position. For ease of exposure, 
we describe each possible position as the set of the sensor 
nodes covered by it. For example, possible positions in Fig. 
[T]can be formulated as Pi = {x 2 ,X 3 }, P 2 = {xi, 0 : 2 ,X 3 }, 
P3 = {a:i,a;2}, P4 = {xi^x^}, P5 = {a;3,X4}, Pq = {xs}. 

In the second step, the RNSC problem searches a smallest 
subset which covers all sensor nodes, from the set of possible 
positions P. In other words, the second step of RNSC prob¬ 
lem reduces to the classic set covering problem for WSNs. 
As the example in Fig. |2] the RNSC problem with X = 
{xi, 0 : 2 , X 3 , a; 4 , ccs} in Fig. [T] is equivalent to searching for 
a smallest subset covering X of P = {Pi, P2, P3, P4, P5, Pq}. 
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Fig. 3. Impact of the location selection for each relay node. 



Definition 3.3 [Network connectivity problem]. Given a 
set of relay nodes Ycdc distributed in the plane, place a 
smallest set of relay nodes Y/vc such that the undirected graph 
G imposed on Y = Yqdc U Ync is connected. 

The network connectivity problem is also known as the 
SMT-MSP problem ll^ - ll^ . The RNSC problem returns a 
set of possible positions that form the backbone of the low-tier 
WSN. As a possible position is defined as the set of locations 
covering the same group of sensor nodes, we should select a 
specific location in each returned possible position to deploy 
one relay node. As shown in Fig.|3 P = {pi,P 2 :P 3 } is the set 
of returned possible positions and the blue line between two 
locations and rj indicates that \\ri — rj\\ < R. For possible 
positions in P, if we choose locations ri, r2 and to deploy 
relay nodes, then additional relay nodes (e.g., and r4) will 
be needed to maintain the connectivity between r2 and r^. In 
contrast, if we choose locations ri, and r4 to deploy relay 
nodes, then the connectivity of ri, and r4 straightforwardly 
holds. Therefore, the number of relay nodes deployed for the 
connectivity of high-tier WSN highly depends on the deployed 
location of relay nodes for covering the low-tier WSN. To 
save the relay nodes deployed for the connectivity of high-tier 
WSN, the deployed locations of possible positions returned 
by the RNSC problem should be optimized. Section V will 
focus on the optimization of deployed location in each possible 
position. 

IV. Approximation Algorithms for RNSC and 
RNDC Problems 

A. The Local Search Approximation Algorithm for RNSC 
Problem 

Definition 4.1 [Neighbor of a possible position]. For any 

two different elements Pi and Pj in P, if at least one sensor 
node within Pi (or Pj) is covered by Pj (or Pi), Pi and Pj 
are mutual neighbors. 

Let N{Pi) denote the set of neighbors of Pi and Ri = 
(Pi U P 2 U ... U Pfc) — Pi denote the sensor nodes covered 
by N{Pi) but not covered by Pi. As shown in Fig. 0] if Pi 
corresponds to area covered by the green circle, then N{Pi) 
represents the area covered by black circles and Ri represents 
the set of the red points. Then, Xi = X — Pi — Ri corresponds 
to the purple points located in blue region in Fig. 0] 


Fig. 4. The LSAA illustration. The sensor nodes are represented by colorful 
points, and the coverage of relay nodes with different possible positions are 
shown by circles with different colors. 


Definition 4.2 [Group]. For any possible position, say Pi, 
the set of sensor nodes covered by the union of Pi and N{Pi) 
is defined as the group of Pi, denoted by group i. 

Definition 4.3 [Bridge]. The possible position covering 
sensor nodes of group i and group j simultaneously is called 
a bridge between group i and group j. 

We denote A* as a Minimum Set Cover (MSC) of the set 
system {X,P). The set consisting of all bridges belonging to 
A* is denoted by bridge a* ■ 

Given a set of sensor nodes X = {xi,X 2 ,Xn}, let group 
i (consisting of the red points and the black points in Fig. 
0]l denote an arbitrary group of X. As shown in Fig. 0] the 
bridges (the blue circles) of group i cover sensor nodes both 
in group i and other groups, i.e., group i is associated with 
other groups by the bridges of group i. A* may include some 
of the bridges, if we know these bridges belonging to A* of 
group i and remove these bridges, we just need to consider the 
possible positions (the black circles and the green circle in Fig. 
01l which only cover group i to cover group i. In other words, 
group i becomes independent from other groups. Therefore, 
as we remove all the bridges (i.e., bridge a-^) belonging to 
A* of each group, all the sensor nodes are divided into 
several independent groups. But this method is based on the 
assumption that the bridgeA* is known, and bridgeA* should 
be built based on A*, whose solution is NP-complete. Hence, 
in the first phase of LSAA, we just group the sensor nodes 
without removing bridge a"- ■ 

The detail of the first phase of the LSAA is described as 
follows; LSAA starts with the multi-stage grouping phase. At 
each stage of the grouping phase, LSAA picks a possible 
position Pi which covers the greatest number of remaining 
sensor nodes. Then, LSAA allocates the sensor nodes covered 
by Pi and N{Pi) as one group. The grouping procedure above 
repeats until all sensor nodes are grouped. 

Since fewer sensor nodes may lead to less searching time 
and favorable performance, we further classify the sensor 
nodes of group i into two kinds, i.e., Ri (the red points in 
Fig. 01 and Pi (the black points in Fig. 01i. Then, we search a 
SC Mi for Ri. Let Wi denote the number of sensor nodes. 
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which are not covered by Mi, in Pi. If Mi fully covers 
Pi, we have Wi = 0, otherwise, we have Wi > 0. The 
second phase of LSAA is to carry out a local search for 
each group. Specifically, LSAA checks the weight Wi of Pi. 
If Wi > 0, LSAA takes {Pi}UMi as the local search result; if 
Wi = 0, LSAA takes Mi as the local search result. The whole 
procedure is repeated until all sensor nodes are covered. The 
LSAA is summarized in Algorithm [T] 


Algorithm 1 Local Search Approximation Algorithm (LSAA). 

Require: 

A set of sensor nodes X = {xi, X 2 , ■. ■ ,Xn}, wireless 
communication radius of sensor node r and relay node R 
{R = 2r) , and a set of possible positions P. 

Ensure: 

A set of relay nodes Y = {yi, y 2 , ..., Um}, which ensures 
that each sensor node is covered by at least one relay node. 
1 : PP = 0; %0 is the empty set. 

2: U = X- 
3: AP = P-, 

4: while {U ^ 0) do 

5: Sort AP in descending order according to the number 

of remaining sensor nodes covered by each possible 
position; 

6 : Pi =Pop(AP); % “Pop” denotes a stack popping 

operation 
7: hnd iV(P,); 

8 : search Mi by a local search approximation algorithm; 

9: calculate wf, 

10 : switch(w;i) 

11 : case Wi > 0: PP = {Pi} U Mi, 

12 : case Wi = 0: PP = Mi, 

13: j = 1; 

14: while (? < \PP\) do 

15: AP = AP- PP; 

16: TmpSens =sensor nodes covered only by PPf, 

17: fc = 1; 

18 : while (fc < \AP\) do 

19: if {{APk fully covers TmpSens) && (APk covers 

more sensor nodes than PPj)) then 
20: PPj = APk, 

21 : end if 

22: k = k + 1; 

23: end while 

24: j = j + 1; 

25: end while 

26: Tmp =the sensor nodes covered by PP; 

27: U = U — Tmp; 

28: AP = AP - PP; 

29: Y = YU PP; 

30: end while 
31: return Y; 


Before performing the performance analysis of LSAA, we 
have to first give out the definitions of some important symbols 
that will be used in the following proof. Let Ci and Ai denote 
an MSC of group i and the subset (covering the sensor nodes 


of group i) of A*, respectively. Let Bi = Ai — bridge a* 
denote the set difference between Ai and bridge a* ■ Obviously, 
\Ai\ > \Ci\ holds. Let \OPAi\ denote the cardinality of 
a feasible SC (returned by the approximation algorithm) of 
group i. Let = |OPAi|/|Ci| denote an approximation ratio 
of the local search algorithm at group i. Let rmax and rmin 
denote the largest approximation ratio and the smallest ap¬ 
proximation ratio, respectively, i.e., rmax = arg max and 

l<2<m 

7’min = arg min r^. Let Di denote the set of possible positions 

l<2<m 

(which cover the sensor nodes in group i) of bridge a* ■ Let 
Dmax denote the set with largest cardinality among P^s, i.e., 
\Dmax\ e arg max \Di\. 

l<2<m 

Lemma 1: LSAA yields an approximation ratio of (1 + 
where e = m\Dmax\/\A*\. 

Proof: According to the definition of Bi, we have 


Ai = BiU Di, 1 < i < m. ( 1 ) 

Then, the following inequality ^ holds 


\Bi U Di\ > \Ci\, I <i<m, (2) 

where (Ell straightforwardly follow \Ai\ > \Ci\. 

m 

As IJ \Di\ = bridge A*, we have 
2=1 




2=1 


Ua 


> \bridgeA-\- 


(3) 


Thus, there exists a non-negative constant C such that 


^ |Pj I = C -f \bridgeA- \ ■ (4) 

2=1 

Then, we have 

m m m 

\B,,\ + \bridgeA-\ + C 

2=1 2=1 2=1 

m 

2=1 
771 

>c.Ei^*i’ 

i=l 

where ” = ” and ” >„ ” follows from ( 01 ) and (El, respectively. 
According to the dehnition of Ai, we can easily get 


\b,\Jd, 


(5) 


'^\Bi\ + \bridgeA*\ = \^*\- ( 6 ) 

2=1 

Plugging (121 into (0 yields 

m 

\A*\>J2\C^\-C 

2=1 

m m 

^Ei^*i-Eiai (7) 

2=1 2=1 

m 

> ^|Ci| -m|P 

max I ■ 

2=1 
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Fig. 5. A case demonstrating the inefficiency of the Weighted Greedy 
algorithm. All the circles are returned by the Weighted Greedy algorithm, 
where the optimal solution only includes the red circles. 


or 

m 

m\Dmax\ ^ 

|y4*| - |A*| ■ 

Equivalently, 

m 

m\D^ax\ I , ^ V 

|A*| n\A*\ - r^ax\A*\ ■ 

Then, the approximation ratio of LSAA is given by 


( 8 ) 


(9) 


E\OPA\ in 
^ ^m|D. 


I A* 


I A* 


“t" 1)^77103; - (1 “f 


( 10 ) 


which completes the proof. ■ 

Remark 1: Inequality (O shows that the approximation 
ratio of LSAA only depends on Tmax and e. Thus the ap¬ 
proximation algorithm yielding small rmax and e is favorable. 

f;iopA,i 

Lemma 2: The ratio of — has a lower bound of 

^min • 


B. The Algorithm for Local Search 

K. Ali et. al. lH] proposed a Weighted Greedy algorithm, 
which selects the possible position with the largest weight. The 
weight of possible position Pi is defined as; weight{Pi) = 
\U n Pi\ — Q!(|Pi| — |17 n Pi\), 0 < a < l/n, where U 
denotes the set of remaining sensor nodes. A case shown in 
Fig. |5] demonstrates the inefficiency of the Weighted Greedy 
algorithm. Since the possible position P 2 has the largest weight 
after Pi is selected, P 2 will be selected next and another two 
possible positions will be selected to realize the full cover to 
the sensor nodes. As a result, the SC returned by the Weighted 
Greedy algorithm includes 4 possible positions (all circles in 
Fig-Ell, whereas the MSC for this case includes only 3 possible 
positions (red circles in Fig. |5]l. To avoid the inefficiency 
of Weighted Greedy algorithm, a Neighbor First Weighted 
Greedy Algorithm (NFWGA) is proposed. The NFWGA does 
not only distinguish the remaining sensor nodes U from the 
already covered sensor nodes, which are denoted by A—C7, but 
also classify U into two kinds of sensor nodes. The first kind 
of U is denoted by A^v, which is covered by the neighbors 
of the deployed relay nodes. And the second kind of U is 
denoted by U — Xm- The weights of sensor nodes of Ajv and 
U — Xn are denoted by a and /3, respectively. In order to 
increase the node degree of sensor nodes, we consider giving 
the sensor nodes in A — C7 a weight 7 , where the node degree 
of one sensor node is defined as the number of relay nodes 
covering this sensor node in this paper. The weight 7 is set 
far less than a and /3 so as not to impact the size of the set 
cover of the remaining sensor nodes. The weight of possible 
position Pi is calculated as 

weight{Pi) = bxa + ex/3 + Cx 7 , (14) 

where <5, e and C denote the number of sensor nodes that 
are covered by A^r, U — Xn and X — U, respectively. The 
proposed NFWGA is summarized in Algorithmic] 


Proof: According to the definition of A*, we have 


|A*|< 


2=1 


- X] 




According to the definition of Ci, we have 


|y4* I < 'Y2 


Then, we have 


IO-PA7I 


> r„ 


( 11 ) 


( 12 ) 


^ |A*, 

2=1 ' ' 

This completes the proof of Lemma |2] 

Remark 2: Inequality (El implies I 
improve the lower bound of the perfori 
adopt local search algorithms with a low approximation ratio. 


Algorithm 2 Neighbor First Weighted Greedy Algorithm 
(NFWGA)._ 

1 : P = 0 ; %0 is the empty set. 

2 : U = X- 
3: A = 0; 

4: AP = P- 
5 : while [U f 0 ) do 

6 : SN= the remaining sensor nodes covered by the neigh¬ 

bor possible positions of T; 

7 : for all Pi G AP do 

weight[Pi) = 5{SN) x Q;-|-e(S'A) x 13 + (^{SN) x 7 ; 

end for 



10 

P = arg max wei 

(13) 

11 

2 = 1 ,... ,772 

u = u-p-. 


12 

Y = AU {P}; 

■ 

13 

AP = AP - {P}; 

y to 

14 

end while 

15 

return A ; 
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Fig. 6. The illustration of the largest coverage of a group i. The lai'gest 
coverage of Pi and its neighbors is denoted by the big circle, which is an 
incircle of a square with side length of 6r. Obviously, the square can be 
divided into 9 small squares with side length of 2r. According to ED ED. 
each squai'e with side length of 2r can be fully covered by 4 circles with radius 
r. This leads to that the square with side length 6r can be fully covered by 
36 circles with radius of r. 


C. The Time Complexity Analysis of LSAA 

Let ti and L (1 < i < m) denote the running time of 
the NFWGA for group i and the running time of grouping i, 
respectively. Then, the overall running time of LSAA (Tj^saa) 
is calculated as follows 

m 

Tlsaa = [tj + ii). (15) 

i=l 

The complexity of sorting AP as shown in step 5 of 
Algorithm [T] is 0(|J^p). Since the number of relay nodes to 
cover each group is no more than 36 as illustrated in Fig. |6] 
the complexity of the inner loop (steps 14-25) of Algorithm 
[U is 0(|Arp|J^|). As the time complexities searching for Pi 
and N{Pi) are 0(1) and O(IA'plJ^I), respectively, the time 
complexity of grouping i is given by 

u = 0(1 J-|") -f 0(|X|2 |J-|) + 0(1) + 0(|X|2 |J-|) 

= 0(|J-|") + 0(|X|"|J-|). 

For the given sensor nodes X = {a;i,a:2, ■.■,Xn}, there are 
at most |X|(|Ar| + 1) possible positions, and thus we have 
|J^| = 0(|Arp), which together with (fTbl l implies that ti = 
0(|X|4). 

In NFWGA, the time complexity of calculating the weight 
of each possible position is 0(|2fp), and the time complex¬ 
ity of finding the possible position with maximal weight is 
0(|J^|). As the iteration number of the loop is less than \X\, 
the time complexity of NFWGA is given by 

f, = 0(|X|3|J-|) = 0(|X|"). (17) 

According to equation (ffSl l. we have 

m 

Tlsaa-E( 0(l^l") +0(1^1')) 

7^1 (18) 
= m(0(|X|^)). 

Since m < |X|, the time complexity of LSAA is finally 
given by 0(|xh. 


D. The Algorithm for RNDC Problem 

To achieve robustness against the failure of deployed relay 
nodes, the RNDC problem is considered and the LSAA for 
Double Cover (LSAADC) is proposed to solve the RNDC 
problem. LSAADC consists of two steps. Specifically, in Step 
1 the LSAA is run for the input sensor nodes, and in step 2 
for the sensor nodes whose node degrees are less than 2 after 
step 1. The detail of LSAADC is shown in Algorithm |3] 


Algorithm 3 LSAA for Double Cover (LSAADC). 

Require: 

A set of sensor nodes X = {a;i,X2 ,... ,Xn}, the trans¬ 
mission distance for sensor node and relay node is r and 
relay node R {R = 2r) , respectively, and a set P of all 
the possible position. 

Ensure: 

A set of relay nodes Y = {yi,y 2 , which makes 

each sensor node is covered by at least two relay node. 

1: RNSC = 0; 

2 : Y = 0; 

3: SN = 0; 

4: RNDC = 0; 

5: RNSC =the relay nodes returned by LSAA(Ai); 

6: SN =the sensor nodes covered by less than 2 relay nodes; 

7: RNDC =the relay nodes returned by LSAA(S'W); 

8: y = RNSC U RNDC; 

9: return Y ; 


Remark 3: The worst case of LSAADC is that the node 
degree of each sensor node is 1 after the first step of LSAADC 
and all the sensor nodes are taken as the input of the second 
step. Therefore, the running time of LSASDC T^saadc is 
upper bounded by 2 Tlsaa and Tlsaadc = 0(|2 f |®). 

V. Network Connectivity Problem 

The proposed LSAA outputs a set of possible positions as 
the result for the given RNSC problem. Flowever, a possible 
position includes numerous locations. Thus, we should select 
an exact location for each returned possible position to deploy 
the relay node. Consequently, different selections will lead 
to distinct topologies of the high-tier WSNs, therefore, the 
quantity of relay nodes deployed to build the high-tier network 
connectivity varies with the selection. Thereby, the selection 
yielding a high-tier topology with fewer deployed relay nodes 
for the network connectivity is desired. 

In Ha ED-E], the intersection of two circles, whose cen¬ 
ters are two sensor nodes and radius is the communication ra¬ 
dius of sensor nodes, is usually selected as the relay deployed 
location. Inspired by the fact that the denser network usually 
needs fewer relay nodes to build the network connectivity, this 
paper proposes a Relay Location Selection Algorithm (RLSA) 
to minimize the Euclidean distance between the deployed relay 
node and the sink node. As a result, the high-tier WSN shrinks 
towards the sink node and the average distance between any 
two nodes generally decreases (Although this conclusion does 
not hold mathematically, it has been verified by extensive 
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Fig. 7. An illustration on the plane partitions out of the union of geometric Fig. 8. The illustration of the Case (1) in Lemma|3] 
discs a and b. 


simulations). Further, the transmission power and the number 
of newly added relay nodes for maintaining the connectivity 
of the high-tier WSN are reduced. 

For the special case that one possible position only covers 
one sensor node, we can easily conclude that the nearest point 
in the disc (sensor’s communication range) to the sink node 
is the intersection between the line segment (which connects 
the sink node and the sensor node) and the circle boundary of 
disc. However, when the special case does not hold, the relay 
location selection will be very complicated. In this section, 
we first consider a simple case that one possible position only 
covers two sensor nodes and then extend the result to the 
general case that one possible position covers more than two 
sensor nodes. 

A. A Simple Case 

Two sensor nodes are located at points Xi and Xj (1 < 
hj < 7i), and \\xi — XjW < 2r. The communication ranges 
of Xi and Xj are represented by discs a and b, respectively. 
The edge circles of two discs intersect at points C and D. 
Let 1 1 and 1 2 be the tangents of a and b at C, respectively. 
Let I3 and I4 be the tangents of a and b at D, respectively. 
As shown in Fig. |7] the plane out of the union of a and b 
is divided into 8 partitions (i.e., partition I — VIII) formed by 
lines /i, I2, I3, Ia, XjC, XiC, XjD, and XiD. The union of arcs 
CNiD and CN 2 D is termed as the intersection arc of a and 
b, where Ni and N 2 are the two intersections of line XiXj and 
the intersection arc. Let S denote the sink node which may 
be located in any partition from partition I to VIII. Next, we 
will search for a point nearest to S on the intersection arc. Let 
the angle between li and XiXj be denoted by a, and ZCxjXi 
be denoted by /3 as shown in Fig. |7] 

The following Lemma [3] and Lemma 0] tell that the nearest 
position to S on the intersection arc only depends on the 
location of S and the size relationship between a and [3. For 
this, four cases will be discussed in the proofs of Lemma [3 
and Lemma 01 

Case (1); the sink node S lies in partition I or V; 

Case (2); the sink node S lies in partition II, IV, VI, or 
VIII (a < /3); 


Case (3): the sink node S lies in partition VII or III; 

Case (4): the sink node S lies in partition II, IV, VI, or 
VIII {a > P). 

Lemma 3: No matter Case (1) or Case (2) holds, the point 
nearest to S on the intersection arc is always the intersection 
of the line segment xjS or XiS and the intersection arc. 

Proof: We first consider Case (1) and assume that S lies 
in the partition I (i.e., the area between tangents li and I 3 , a s 
shown in Fig. [Hi. For an arbitrary given point T on CN 2 D 
e xcep t C and D, we know that the line segment ST intersects 
CNiD at one point denoted as T*. Obviously, T* is closer 
to S than T. Therefore, the nearest point to S must be on 
CNiD and in the rest of t his p roof we will only search for 
the nearest point to S on CNiD. 

Since the partition I is also located between the hal f lines 
XjC and xjD, the line segment xjS intersects CNiD. Let 
H be the intersection of XjS and CNiD. Next we will show 
that H is the point nearest to S on CNiD. 

Let F be an arbitrary point on CNiD except for H, and I 3 
be the tangent of b passing through H. Evidently, I 3 _L XjS. 
As S and F are located at the two sides of I 3 , respectively, 
i.e., ZFHS > 90°, according to the law of cosines, we know 
that FS is the longest side of AFFIS, i.e., FIS < FS. In 
other words, H is the nearest point to S on the intersection 
arc. Because of the symmetry of discs, the case when S lies 
in the partition V can be similarly proved and thus skipped. 
So far, we have proved the first half of Lemma [3] 

Then we consider Case (2) and assume that S lies in the 
partition VIII (i.e., the area between the half line XjC and 
li , as shown in Fig. |9]l- Let G denote an arbitrary point on 
CN 2 D. Since a < ^ and S lies in the partition VIII, we 
can see that S lies a bove the line li, which implies that SG 
may i ntersect with CN 2 D directly without intersecting with 
CNiD. Therefore, we cannot prove the case (2) by straightly 
following the proof method of case (1). 

Applying the law of cosines to ZSxiG and ZSxiC, we 
have that 


SC^ = -\- XiS'^ — 2 rxiS cos ZSxiC 

SG^ = r'^ + XiS'^ — 2 rxiS cos ZSxiG. 


( 19 ) 
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Fig. 9. The illustration of the Case (2) in Lemma[3] 



Fig. 10. The illustration of the Case (3) in Lemma |4] 


Since the cos(a;) decreases in x when 0 < a; < tt, we conclude 
that S C < SG, which implies that C is the nearest point to 
S on Cl^. 

Since a < /3 and S is located under the line Xj S, obviously, 
XjS and CNiD intersect at a point E. By the similar method 
i n Ca se (1), we can prove that E is the nearest point to S on 
CNiD, which implies that E is closer to S than C and thus 
is the nearest point to S on the intersection arc. Because of 
the symmetry of partitions, the case when S is located in the 
partition II, IV or VI (a < /?) can be similarly proved and 
thus skipped. This completes the proof of Lemma |3] ■ 

Lemma 4: No matter Case (3) or Case (4) holds, the point 
nearest to S on the intersection arc is always the point C or 
D. 

Proof: First of all, we consider the Case (3). Again, 
because of the symmetry of discs, we only consider the case 
that S is located in the partition VII (i.e., the area between 
the half lines XjC and XiC, as shown in Fig. [Toll . 

Let G denote an arbitrary point on the minor arc CNi 
except for C. Let Iq be a parallel line (passing through G) 
of line XiXj, i.e., Iq || XiXj. We then draw an orthogonal line 
GE of le at G (i.e., GE _L Iq), where E is the intersection of 
GE and DNi. As S and E are located at the two sides of Iq, 
i.e., ZSGE > 90°, according to the law of cosines, we know 
that SE is the longest side of ASGE, which implies that the 
nearest point to S can never on the minor arc IDNi. Similarly, 
we can skip the minor arc ]DN 2 and only search the nearest 
point to S on minor arcs CNi and CN 2 . 

Next we search the nearest point to S from the minor arc 



Fig. 11. The illustration of the Case (4) in Lemma |4] 


CNi. Applying the law of cosines to ZxjCS and ZxjGS 
yields 

SC'^ = + XjS'^ — 2rXjS cos ZSxjC 

SG'^ = + XjS'^ — 2rXjS cos ZSxjG. 

Since the cos(x) decreases in x when 0 < x < tt, we can 
conclude that SC < SG since ZSxjC < ZSxjG, which 
implies that C is nearest point to S on the minor arc CNi. 
Similarly, based on the law of cosines, we can achieve the 
conclusion that C is the nearest point to S on the minor arc 
CN 2 . Thus, we have proved that C is the nearest point to 
S on the intersection arc when S is located in the partition 
VII. For the symmetry of discs, we can conclude that D is 
the nearest point to S on the intersection arc when S lies at 
the partition III. 

Then, we consider the Case (4). We take the case when S is 
located in the partition VIII as shown in Fig.[TT]as an example. 
Since a> jS, partition VIII is the area between the half lines 
XjC and XiC, which implies that neither XjS can intersect 
with CNiD nor XiS can intersect with CN 2 D. Let G be an 
arbitrary point except C on the intersection arc. SC < SG 
directly follows from the same proof for Case (3). Therefore, 
C is the nearest point to S on the intersection arc when S 
lies in partition VIII. Finally, for the symmetry of discs, the 
case when S is located in the partition II, IV, VI and a > /?, 
can be proved through the same method. This completes the 
proof of Lemma |4] ■ 

Combining Lemma 0 and Lemma |4] we conclude that the 
nearest point to S on the int ersec tion arc is the inters ectio n 
point if either XiS intersects CN 2 D or XjS intersects CNiD, 
and otherwise, the nearest point to S on the intersection arc 
is C or D. The above analysis is summarized in Algorithmic 

B. The General Case 

This subsection extends the result of the simple case to the 
general case. The proposed RLSA is a multi-stage algorithm, 
with each stage consisting of three steps. In the first step 
of each stage, we randomly select one possible position P 
and search the candidate locations (including the intersections 
and nearest point to S) on the intersection arc of each pair 
of sensors covered by P. Then, we select the locations fully 
covering P among those candidate locations returned by the 
first step. Finally, we find a nearest location to S by extensive 
search among these locations selected in the second step. The 
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Algorithm 4 Nearest Point to the Sink (NPS). 

Require: 

A pair of discs a and b with radius r, centered at xi and 
X 2 , respectively, ||xi — a; 2 || < 2r, the position of the sink 
node S. 

Ensure: 

A point p that is nearest to S on the intersection arc of 
a and b. 

1 : p = 0; 

2: cnt = 1; 

3: find the intersections, which are denoted by C and D, of 
xi and X 2 ', %we denote the edge circles of a and bhy xi 
and X 2 , respectively. 

4: while {cnt < 3) do 

5: p =the intersection between Sxcnt and Xcnt', 

6: if {p lies on the intersection arc) then 

7: break; 

8 : end if 

9: p = 0; 

10: cnt = cnt + 1; 

11 : end while 

12 : if (p == 0 ) then 

13: p = the point, which is nearest to S among points C 

and D; 

14: end if 
15: return p\ 




• sensor node 


Fig. 12. An example of the general case (four discs). 


procedure above repeats until all returned possible positions 
are deployed. The RLSA is described in Algorithm |5] 

Fig. [12] reloads the example in Fig. [T] and uses it to 
illustrate how the RLSA works. It is evident form Fig. [T2| 
that {FbjFsjFe} is a possible position set (i.e., SC) that 
fully covers the sensor nodes X = {aii, 0 : 2 ,aia, X 4 , 0 : 5 }. First, 
we determine the deployed location of P 2 = {xi,X 2 ,X 3 }. 
According to the RLSA, in the first step we search the 
intersections and nearest points for each pair of xi, X 2 , X 3 , and 
these candidate locations are denoted by si, S 2 , ■■■, sg. In the 
second step, we find that only { 53 , 54 , 35 , 56 } (denoted by the 
blue squares and red triangles) can fully cover {xi,X 2 ,X 3 }. 
Finally, we search the nearest point 53 (denoted by the red 


triangle) to S from (53, 54, S5, 55} by exhaustive search. As P5 
has only two candidate locations (siq and sn), 5io (denoted by 
the red square) is finally selected. Since Pq has no neighbors, 
the intersection 5ii (denoted by the red triangle) of the line 
segment x^S and the circle centered at X 5 is selected as the 
deployed location of Pq. So far, we have selected locations 
53, 510 and sii as the relay deployed locations of P2, P5 and 
Pq, respectively. 


Algorithm 5 Relay Location Selection Algorithm (RLSA). 

Require: 

A set of sensor nodes X = {xi,X 2 , ■ ■ ■ ,Xn}, a set of 
possible positions P = {pi,p 2 , ■■■,Pm}, the location of 
the sink node S. 

Ensure: 

A set of locations L = {li,l 2 , to deploy relay 

nodes. 

1: L = 0; 

2 : while (P 7 ^ 0 ) do 
3: TmpP =Pop(P); 

4: TmpX =sensor nodes covered by TmpP of X-, 

5: copyX = TmpX', 

6: while {TmpX ^ 0) do 

7: Tmpl =Pop{TmpX); 

8: neigh =neighbor sensor nodes of Tmpl in TmpX', 

9: if {neigh = 0) then 

10: PPTmp=the nearest point to S on the disc cen¬ 

tered at Tmpl', 

11 : end if 

12: while {neigh 7 ^ 0) do 

13: Tmp2 =Pop{neigh)', 

14: Tmp3 =NPS{Tmpl, Tmp2, S); 

15: TmpA =intersections of Tmpl and Tmp2', 

16: TmpS = Tmpi U TmpA', 

17: PPTmp = TrnpH U PPTmp', 

18 : end while 

19: end while 

20: TmpH =points (of PPTmp) fully covering copyX', 

21: TmpA =point (of Tmpi) nearest to S'; 

22: L = LLIT mpA', 

23: end while 
24: return L; 


In RLSA, the time complexity of the inner loop between 
step 12-18 is 0(1X1) and thus the complexity of the double 
loop between step 6-19 is 0(|Xp). Since the iteration number 
of the main loop is |P|, the time complexity of RLSA is finally 
given by 0(|P||X|2). 

In this paper, we combine the RLSA with the MST heuristic 
proposed in EH to build the network connectivity. Specifi¬ 
cally, when the RNSC problem or RNDC problem is solved, 
we apply the RLSA to search an optimal deployed location 
for each possible position returned by the LSAA or LSAADC. 
Then the MST heuristic is employed to build the connectivity 
for the high-tier network by adding relay nodes if needed. 
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Number of sensor nodes n 
Fig. 13. The influence of a and 0. 


VI. Performance Evaluation 

In the simulation, sensor nodes are randomly deployed in 
a square area with the size 100 x 100 m^. The number of 
deployed sensor nodes n ranges from 10 to 100. Communica¬ 
tion radius for sensor node and relay node are r = 10m and 
R = 20m, respectively. The comparison results are averaged 
over 100 simulations. The simulation platform is developed 
by the C-H-, and simulations are executed on a 3.40-GHz 
Windows Work Station with 16 GB memory. 

A. Parameters Comparison for NFWGA 

Since the NFWGA is taken as the local search algorithm of 
LSAA, we first determine the value of parameters of NFWGA. 
In NFWGA the values of a and /3 determine the size of set 
cover, and 7 is set quite small compared with a and /3 so as 
not to influence the size of set cover. To be Specific, the exact 
value of 7 is set 0.01. For practical use, we obtain the value 
of (a, P) by simulation. Fig, flj] compares the performances of 
NFWGA under different values of (a,/3). Comparison results 
show that the {a, /3) = (5,1) yields the fewest deployed relay 
nodes. Without loss of generality, we set {a, P) = (5,1) in 
the following simulations. 

B. Performance Comparison for LSAA 

To demonstrate the efficiency of FSAA, the algorithms 
based on the Shift Strategy, the Grid Strategy and the Set- 
Covering Strategy are taken for comparison. For simplicity, 
in the rest of this paper, we denote the algorithm proposed 
in csi, ini and ED as Shift, Grid and Weighted Greedy, 
respectively. 

Fig. [ 14 ] compares the FSAA with ifTsll . ifTTl and IJD. 
The local search phase of FSAA adopts the NFWGA. It is 
clearly shown that FSAA returns the least number of relay 
nodes among all four algorithms. We observe that a large 
number of relay nodes could be saved due to FSAA in 
comparison to ifTsl . ifTTl and ll 2 D . in which the largest saving 
(8/24 Ri 33.3%) occurs when n = 40 and the Grid algorithm 
is adopted. 



Number of sensor nodes n 

Fig. 14. Comparisons on the number of deployed relay nodes 

between FSAA and existing works. 



1.5-|| Weighted Greedy [21] 

I ^Blsaa _ I 

mi 

20 40 60 80 100 

Number of sensor nodes n 

Fig. 15. Comparison on the average node degree of sensor nodes between 
LSAA and existing works 


Fig. [15] shows that the Grid algorithm is the most robust 
covering algorithm among the four algorithms. However, the 
poor performance in the number of deployed relay nodes limits 
the use of the Grid algorithm in practice. Although fewer relay 
nodes are deployed by FSAA, the average node degree of 
FSAA is generally larger than the Weighted Greedy algorithm, 
and comparable to the Shift algorithm as shown in Fig. [T5| 
which can be explained by the introduction of weight 7 in 
NFWGA. To be specific, when multiple possible positions 
have the same weight 5xa-lex/3on remaining sensor nodes, 
the possible position covering more covered sensor nodes will 
be selected. 

Fig. [T6| shows that the running times of four compared 
algorithms increase as n increases. FSAA has a moderate 
running time among the four algorithms. The largest running 
time of FSAA (happening when n = 100) is around 275ms. 
Therefore, the time overhead of FSAA is affordable for 
practical use. 

C. Performance Comparison for LSAADC 

The algorithm proposed in lITSl to the RNDC problem 
is denoted by ShiftDC in this paper. Fig. [17] shows that 
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Fig. 16. Comparison on the running time between LSAA and existing works. 


60 n 



20 40 60 80 100 


Number of sensor nodes n 

Fig. 17. Comparison on the number of deployed relay nodes between 
LSAADC and ShiftDC. 


LSAADC deploys fewer relay nodes for double cover than 
the ShiftDC. Meanwhile, in comparison to ShiftDC a slightly 
higher average node degree is achieved by LSAADC as shown 
in Fig. [18] 

In addition, the comparison on running time between 
LSAADC and ShiftDC is shown in Fig. [T9| Similar to the 
comparison between LSAA and the Shift algorithm in Fig. 
[TbI ShiftDC also outperforms LSAADC and the performance 
gap between ShiftDC and LSAADC becomes more visible in 
comparison to Fig. [Tb] 

D. Performance Comparison for RLSA 

In order to examine the effectiveness of RLSA, we first 
show that the high-tier network will become dense when 
we perform RLSA to select relay deployed locations. Fig. 
l20l compares RLSA with the Intersection Location Selection 
Method (ILSM) 1(151 lISTl and a Random Location Selection 
Method (RLSM) in terms of the average distance between each 
pair of deployed relay nodes. We can see from Fig. |20| that the 
RLSA returns the smallest average distance among all three 
algorithms and RLSM returns the largest average distance. 

Fig EHthen presents the comparison on the average number 
of deployed relay nodes between RLSA-based, ILSM-based 



20 40 60 80 100 

Number of sensor nodes n 


Fig. 18. Comparison on the average node degree of sensor nodes between 
LSAADC and ShiftDC. 



Fig. 19. Comparison on the running time between LSAADC and ShiftDC. 


and RLSM-based connectivity algorithms. Notice that, for 
fair comparison, all of the three compared algorithms adopt 
an MST-based algorithm-Feasible single-tiered Relay Node 
Placement (FltRNP) 123 to solve the network connectivity 
problem. As RLSA selects the nearest deployed position for 
each relay node, the average distance between relay nodes 
and the sink node of RLSA will be the smallest. In other 
words, because of RLSA the high-tier network will shrink 
towards the sink node. This explains that RLSA-based network 
connectivity algorithm needs fewer relay nodes for network 
connectivity than the other two algorithms. It is shown in Fig. 
im that the RLSA-based connectivity algorithm at most saves 
13.3% deployed relay nodes comparing with other algorithms 
when n = 60. 

Fig. |22| shows the comparison on the running time of RLSA, 
ILSM and RLSM. It is clearly shown that as n increases, 
the running time of RLSA grows much faster than those of 
ILSM and RLSM, which can be explained by the fact that the 
inner loop (Step b-19) of Algorithm 5 (RLSA) contributes an 
additional time complexity O(n^). 

Finally, in order to evaluate the overall performance of this 
work (denoted by LSAA-RLSA, which consists of LSAA 
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Number of sensor nodes n 


Fig. 20. Comparison on the average distance between each pair of deployed 
relay nodes. 



Fig. 23. Comparison on the average number of deployed relay nodes between 
LSAA-RLSA, Shift-ILSM, Grid-ILSM 

and Weighted Greedy-ILSM. 




Fig. 24. Compai'ison on the running time between LSAA-RLSA, Shift-ILSM, 
Grid-ILSM and Weighted Greedy-ILSM. 


Fig. 21. Comparison on the average number of deployed relay nodes between 
RLSA-based, ILSM-based and RLSM-based connectivity algorithms. 



Fig. 22. Comparison on the running time between RLSA, ILSM and RLSM. 


and the RLSA-based network connectivity algorithm), the 
algorithm presented in ifTSl . ifTTl and IJl], which are denoted 
by Shift-ILSM, Grid-ILSM and Weighted Greedy-ILSM, re¬ 
spectively, are taken for comparison. All of the three compared 
algorithms employ FltRNP EH to solve the network con¬ 
nectivity problem. Fig. |23 and Fig. |24] show the number of 
deployed relay nodes and running time of these algorithms, 
respectively. To sum up, LSAA-RLSA is the most time- 
consuming but at the same time also the most cost-effective 
algorithm. As the largest running time of LSAA-RLSA in 
simulation is 550ms (which is still affordable for practical 
use), the proposed LSAA-RLSA is the most favorable one 
among the four algorithms because of its superiority in the 
cost of relay deployment. 

VII. Conclusion 

In this paper, we have investigated the relay node place¬ 
ment problem in WSNs. Firstly, by modeling the relay node 
placement problem as a GDC problem, we proposed a novel 
algorithm-LSAA to solve it. The novelty of LSAA lies in the 
separation of the grouping phase and the local search phase. 

















































































14 


The LSAA yields an approximation ratio of (1 + e)ri„ax, 
and the time complexity of LSAA is proved to be 0(|X|®), 
where rmax denotes the maximal local approximation ratio. 
The number of deployed relay nodes of LSAA is compared 
with the existing works through simulations and the results 
show that at most one third of relay nodes could be saved 
due to this new approach of LSAA. Secondly, in order to 
improve the network robustness against connection failure, the 
RNDC problem was solved to enhance the node degree of the 
network. The efficiency of the solution algorithm-LSAADC to 
the RNDC problem is demonstrated through simulations. The 
simulation results show that the LSAADC deploys fewer relay 
nodes while achieving a higher node degree in comparison to 
existing works. Thirdly, a unique RLSA-based NC algorithm 
was proposed to search a nearest deployed location to the 
sink node for each relay node. Then the MST heuristic is 
employed to build the connectivity for the high-tier WSN. 
Again, the number of newly added relay nodes for building 
the connectivity of the high-tier WSN is significantly saved 
in simulations. In summary, based on comparisons with the 
published works and the presented results, this work outlines 
a new path for relay node placement for the new generations 
of WSN applications. 
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